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Abstract 

Assuming that an unknown mechanism (e.g., gas turbulence) removes most 
of the subnebula gas disk in a timescale shorter than that for satellite formation, 
we develop a model for the formation of regular (and possibly at least some 
of the irregular) satellites around giant planets in a gas-poor environment. In 
this model, which follows along the lines of the work of Safronov et al. (1986), 
heliocentric planetesimals collide within the planet's Hill sphere and generate 
a circumplanetary disk of prograde and retrograde satellitesimals extending as 
far out as ~ Rh/%- At first, the net angular momentum of this proto-satellite 
swarm is small, and collisions among satellitesimals leads to loss of mass from 
the outer disk, and delivers mass to the inner disk (where regular satellites 
form) in a timescale < 10 5 years. This mass loss may be offset by continued 
collisional capture of sufficiently small < 1 km interlopers resulting from the 
disruption of planetesimals in the feeding zone of the giant planet. As the 
planet's feeding zone is cleared in a timescale < 10 5 years, enough angular 
momentum may be delivered to the proto-satellite swarm to account for the 
angular momentum of the regular satellites of Jupiter and Saturn. This feeding 
timescale is also roughly consistent with the independent constraint that the 
Galilean satellites formed in a timescale of 10 5 — 10 6 years, which may be long 
enough to accomodate Callisto's partially differentiated state (Anderson et al. 
1998; 2001). In turn, this formation timescale can be used to provide plausible 
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constraints on the surface density of solids in the satellitesimal disk (excluding 
satellite embryos ~ 1 g cm" 2 for satellitesimals of size ~ 1 km), which yields a 
total disk mass smaller than the mass of the regular satellites, and means that 
the satellites must form in several ~ 10 collisional cycles. However, much more 
work will need to be conducted concerning the collisional evolution both of the 
circumplanetary satellitesimals and of the heliocentric planetesimals following 
giant planet formation before one can assess the significance of this agreement. 
Furthermore, for enough mass to be delivered to form the regular satellites 
in the required timescale one may need to rely on (unproven) mechanisms to 
replenish the feeding zone of the giant planet. We compare this model to 
the solids-enhanced minimum mass (SEMM) model of Mosqueira and Estrada 
(2003a,b), and discuss its main consequences for Cassini observations of the 
Saturnian satellite system. 
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1 Introduction 



The similarities between planetary and satellite systems led to early suggestions that 
any theory for planetary formation should incorporate regular satellite formation as 
a natural byproduct (Alfven 1971; De et al. 1977). More recently the propensity has 
been to argue against a unifying, generally applicable formation model for planets 
and satellites {e.g., Stevenson et al. 1986). Yet, whether or not regular satellites and 
planets formed in an analogous manner, it is often remarked that the scaled bulk 
characteristics of angular momentum and mass are shared by the regular satellites 
of Jupiter, Saturn and Uranus {e.g., Pollack et al. 1991). The first to suggest that 
regular satellites formed out of a collisionally-captured gravitationally-bound swarm 
of circumplanetary satellitesimals was O. Yu Schmidt (1957). This idea was subse- 
quently explored in a number of publications (Safronov and Ruskol (1977); Ruskol 
1981, 1982; Safronov et al. 1986). However, it is fair to say that while these works pro- 
vide valuable insight into the physical mechanisms possibly involved in the formation 
of regular satellites of giant planets out of an extended circumplanetary satellitesimal 
swarm, they fall short of advancing a specific model that can account for the observed 
properties of the satellite systems. In this contribution, we revisit this approach with 
the aim of providing a testable model as well as a basis for comparison with our pub- 
lished solids-enhanced minimum mass (SEMM) satellite formation model (Mosqueira 
and Estrada 2003a,b; hereafter, MEa,b). 
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The primary difference between the SEMM model and the present gas-poor plan- 
etesimal capture (GPPC) model is in the treatment of turbulence. It is commonplace 
in the planetary formation literature to assume that Keplerian disks are turbulent 
and to model the resulting anomalous angular momentum transport in terms of an 
a ~ 1CT 3 parameter 1 (Shakura and Sunyaev 1973) regardless of disk conditions or 
distance from the central object. While this assumption may or may not (given the 
scatter in the data and the existence of alternative explanations) be justified in the 
case of the solar nebula by the observation that disks around young protostars accrete 
at low rates ~ 1O~ 8 M yr _1 (Hartmann et al. 1998), there is no reason to suppose 
that similar a values need also apply to circumplanetary nebulae at the time of satel- 
lite formation. Indeed, even those turbulent mechanisms (such as MHD [Balbus and 
Hawley 1998] and baroclinic instabilities [Li et al. 2000; Klahr and Bodenheimer 
2003]) that may operate at least at some locations and times in the case of the 
planetary disk may fail in the dense (but not massive), cold and mostly isothermal 
circumplanetary disk (MEa). 2 

It is often suggested that the disk's radial shear may drive turbulence for suffi- 
*It should be kept in mind, however, that some workers favor the /3 parameter (i.e., Richard and 

Zahn 1999) 

2 One would expect subnebula turbulence while there is significant Roche- lobe gas inflow, but it 
is difficult to see how this source of turbulence can be relied upon to remove the circumplanetary 
gas disk. 
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ciently high Reynolds number (Zahn 1991; Dubrulle 1993; Richard and Zahn 1999). If 
true, such a mechanism would apply quite generally to accretion disks (both planetary 
and satellite disks would be affected); however, so far numerical simulations (Hawley 
et al. 1999; Mosqueira et al. 2003) have failed to show the presence of a non-linear 
shear instability, apparently due to the stabilizing Coriolis terms in flows with posi- 
tive radial gradients in specific angular momentum (cf. Rayleigh stable flows), even 
in the presence of both a radial and a vertical shear (Rudiger et al. 2002). In fact, 
it has been claimed that a carefully constructed Taylor-Couette experiment shows 
stability in the case of positive radial gradients in specific angular momentum even 
for high Reynolds number (Schultz-Grunow 1959) 3 . The difficulty stems not only 
from the degree of turbulence, but also whether such turbulence, if present, may be 
characterized by an a parameter. 

Given the lack of observational or theoretical quantitative constraints for the gas 

inflow rate through the gap or planetary accretion rate in the final stages of giant 

planet formation (Lubow et al. 1999), characterization of the subnebula evolution at 

the time of satellite formation is clearly beyond reach at the present time. One may 

tackle this hurdle either by assuming that turbulence decays so that Keplerian disks 
3 Though it is unclear whether this would remain true for even higher Reynolds number. Ex- 
perimentally, the onset of turbulence in pipe flows cannot be described simply by linear stability 
analysis. 
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must pass through quiescent phases (which enables planetesimal and satellitesimal 
formation), or by assuming that satellites begin to form once the subnebula has been 
mostly depleted of gas (which facilitates their survival) presumably due to sustained 
turbulence of unknown origin (possibly hydrodynamic shear turbulence). In the for- 
mer case, satellite survival hinges on gap-opening (by satellites with mass ratio to 
the primary fj, ~ 1CT 4 ), and this mechanism itself sets the value for the gas surface 
density of the satellite disk (MEb). In the latter case, which is the focus of the present 
paper, the gas surface density is low but unspecified, though the presence of some 
gas may help to explain the observations. Here satellite formation is dealt with in a 
manner somewhat analogous to the consensus treatment for the terrestrial planets. 

The second major difference between the two models concerns the mechanism con- 
trolling the timescale for satellite formation. MEa,b includes satellitesimal migration 
due to gas drag and tidal torques, and forms satellites by a combination of Safronov- 
type binary accretion and drift augmented accretion in an extended, two-component, 
largely quiescent subnebula. This leads to a model in which satellitesimals (~ 100 
km) form quickly (~ 1000 years) but full grown satellites (~ 1000 km) take signif- 
icantly longer to form (10 4 — 10 6 years depending on their location; Mosqueira et 
al. 2001), in a time set by the gas drag time of satellitesimals in a solids-enhanced 

minimum mass subnebula 4 (enhanced by a factor ~ 10, which leads to a gas surface 
4 Gas flow through the giant planet's gap (Lubow et al. 1999) does not necessarily imply that 
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density ~ 10 4 g cm -2 ). Furthermore, this formation timescale is consistent with the 
Type I migration timescale of full grown satellites in such a disk provided 3-D effects 
are considered (which leads to a factor ~ 10 slower migration; e.g., Bate et al. 2003), 
and with the criterion for gap-opening in an inviscid disk (with aspect ratio ~ 0.1; 
Rafikov 2002). 

In contrast, in this paper the formation timescale for the regular satellites is set 
by the timescale over which planetesimal fragments are collisionally captured by the 
circumplanetary satellitesimal swarm, which is ultimately tied to the timescale for a 
giant planet to clear its feeding zone (Mosqueira et al. 2000). Given that the gas- 
surface density in this model is unspecified but is presumably quite small (< 100 g 
cm -2 ), the migration of satellites due to either gas drag or tidal torques does not pose 
a threat to their survival; yet, one may still use the gas disk to circularize the orbits 
of the regular satellites, or to sweep away collisional debris from regions of the disk 
that are observed to be either empty or significantly mass-depleted (such is the case 
with most of the Saturnian satellite system). 

The third major difference we draw attention to in this introduction (see also 

the comparison table in Appendix C as well as MEa,b) concerns the interpretation 
the time for disk formation exceeds that of satellite formation. For instance, one possible scenario 

is for Jupiter and Saturn to clear the disk in a timescale < 10 5 years in a weakly turbulent nebula 

(Brydcn et al. 1999, 2000). 
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of the compositional gradient of the Galilean satellites. Two main explanations have 
been suggested for this observation. A number of publications ascribe the high-silicate 
fraction of Europa (Io might have lost its volatiles due to the Laplace resonance alone) 
compared to Ganymede and Callisto to the subnebula temperature gradient due to the 
giant planet's luminosity at the time of satellite formation following envelope collapse 
(e.g., Pollack et al. 1976; Lunine and Stevenson 1982). On the other hand, Shoemaker 
(1984) argued that this compositional gradient might be due to the increase of impact 
velocities and impactor flux of Roche-lobe interlopers deep in the planetary potential 
well, leading to preferential volatile depletion in the case of Europa and Io. That is, 
in this view all the Galilean satellites started out ice-rich, but some lost more volatiles 
than others. 

Several more recent papers attempt to explain the volatile composition of the 
Galilean satellites based on turbulent models of the Jovian subnebula. Makalkin et 
al. (1999) identify the difficulties involved in explaining the Galilean satellite com- 
position gradient using turbulent subnebula models, but do not attempt to produce 
a specific model of satellite formation. Mousis et al. (2002) and Mousis and Gau- 
tier (2004) employ an evolutionary turbulent subnebula and deliver solids to form 
satellites by "embedding" solar nebula planetesimals that incorporate clathrate hy- 
drates. However, how heliocentric planetesimals become embedded in the subnebula 
is never addressed, and the mechanism for trapping volatiles is uncertain (Lunine et 
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al. 2000). On the basis of this clathration theory, Hersant et al. (2004) have repro- 
duced the abundances of C, N and S in Saturn available in 2004, and predicted that 
only Xe is substantially enriched in Titan's atmosphere. However, the CIRS instru- 
ment aboard Cassini (Flasar et al. 2005a) has found a C/H ratio about twice higher 
than predicted by Hersant et al. (2004), and no noble gases other than radiogenically 
produced 40 Ar (and trace amounts of 36 Ar) were detected in Titan's atmosphere by 
the GCMS aboard Huygens (Niemann et al. 2005). Hersant et al. (2005) explain 
these discrepancies in terms of the uncertainty in the C/H ratio and argue that Xe 
may remain trapped in Titan's interior. 

Canup and Ward (2002) posit a gas-starved steady-state subnebula such that a 
balance is attained between gas inflow through the giant planet's gap and turbulent 
removal from the circumplanetary disk. In this model a combination of viscous heat- 
ing and heating from Jupiter's luminosity prevent condensation of ices close to the 
planet, which explains the Galilean satellite compositional gradient. However, these 
authors assume that the gas inflow forms a compact disk around planet, which then 
viscously spreads. Yet, the gas inflow through the giant planet's gap would be en- 
dowed with enough specific angular momentum to form an extended circumplanetary 
gas disk (MEa; note that this paper considers gas flow through the gap in connection 
with Callisto), which would result in satellites forming in very cold regions of the disk 
with their full complement of ices. 
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Finally, Alibert et al. (2005) forms satellites in a disk essentially like that of MEa, 

except that they retain disk evolution due to weak turbulence instead of allowing for 

turbulence decay. Like MEa, they employ an inner disk out to ~ 30Rj and an 

outer disk out to ~ Rh/5, and similar inner disk but larger outer disk (note that 

this transition is not as pronounced as in MEa) gas surface densities. Moreover, 

their assumed satellite formation timescale (10 4 — 10 6 yrs depending on location; 

Mosqueira et al. 2001; MEa) is controlled by disk dynamics (as supposed to the gas 

inflow timescale as in Canup and Ward 2002), which means that satellite observations 

would receive very similar explanations to those outlined in MEa. However, unlike in 

MEa,b these initial conditions are not justified. Moreover, too much mass is located in 

the cold, outer regions of the disk, so that one would expect satellites to form far from 

the planet (where they are not found). On the other hand, if such satellites migrated 

large distances it would then become difficult to explain how rocky satellites close to 

the planet managed to survive. Also, these authors base their temperature profiles 

on the treatment of Papaloizou and Terquem (1999), which does not include heating 

by the central object. Thus, Alibert et al. (2005) attempt to explain the Galilean 

satellite compositional gradient using viscous heating alone. As a result, their model 

yields implausibly low temperatures both for the nebula and the subnebula as gas 

dissipation takes place 5 (c.f. their Fig. 6 has a temperature of T ~ 30 K at 10Rj 
5 This assumes shielding by the inner disk. Whether this is feasible remains to be shown. At any 
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and T ~ 10 K at 150Rj for t — 1 Myr). Furthermore, their model relies on a much 
different Type I migration attenuation factor // for the nebula than for the subnebula. 

In the context of our decaying-turbulence SEMM model, the Galilean satellite 

compositional gradient is used to set a subnebula temperature at the location of 

(gap-opening) Ganymede of T ~ 250 K (or aspect ratio of ~ 0.1 as stated above) 

controlled by the giant planet's luminosity. 6 This may imply a planetary luminosity 

between 10~ 4 — 1O~ 5 L for a planetary radius of 1.5 — 2.0Rj (where Rj is Jupiter's 

radius) following envelope collapse (work in progress). Here we simply note that 

deep in the planetary potential well one might expect a stochastic compositional 

component due to high speed ~ 10 km s -1 impacts with large ~ 100 km Roche-lobe 

interlopers (but see Mosqueira and Estrada 2005 for the case of Iapetus, which may 

be exempt from this argument due to its distance from the planet, size and isolation). 

Such hyperbolic collisions might conceivably remove volatiles from the mantle of a 

differentiated satellite and place them on neighboring satellites (possibly analogous 

to the impact that may have stripped Mercury's mantle; Benz et at. 1988). This 

might explain the volatile depletion in Europa relative to its outer neighbours (given 

its density [ Jacobson et al. 2005] , a similar argument may apply to Enceladus in the 

rate, this might again argue against a compositional gradient based on subnebula temperature. 
6 It is important to stress that because this is a decaying turbulence model, the disk's temperature 

is determined by the planet's luminosity and not by viscous heating. 
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Saturnian system; alternatively, a different fractionation mechanism may apply to 
both) and obviate the need for a temperature of formation constraint. In this paper 
we adopt this point of view. In this connection, it is worth noting that recent work 
suggests that close-in Amalthea may have incorporated water-ice (Anderson et al. 
2005), which is a definite possibility in the framework of the present model. 

Finally, it is useful to compare the method of delivery of solids in our two mod- 
els. We expect that the regular satellites form at the tail-end of planetary accretion 
so that, at least in the core accretion model for the formation of the giant planet, 
almost all of the solids mass in the solar nebula starts out in the form of ~ 10 km 
planetesimals (e.g., Weidenschilling 1997). However, following giant planet formation 
fragmentation may cascade a significant fraction of the mass in the feeding zone of 
the giant planets back to smaller sizes (Stern and Weissman 2001; Charnoz and Mor- 
bidelli 2003). Planetesimals in the meter to kilometer mass range may deliver mass 
to the regular satellites of the giant planets either by ablating in the giant planet 
envelope prior to envelope collapse, or by ablating through the circumplanetary gas 
disk (MEa; Mosqueira and Estrada 2005). In the context of this paper, planetesimals 
in this size range are both easier to capture into the circumplanetary disk, and may 
also be used to replenish the feeding zone of the giant planet by gas drag. 

In Section 1.1 we summarize the main points of the GPPC model. The creation of 
the circumplanetary satellitesimal disk is discussed in more detail starting in Section 

12 



2. Readers who are more interested in the results may wish to skip to Section 2.3. 
In Section 2.4, we present results of the angular momentum calculation. Finally, in 
Section 3 we present discussion of the results and our conclusions. Several appendices 
are provided at the end of the paper for reference. 

1.1 Model Outline 

In this section, we summarize the GPPC model. In the spirit of Safronov et al. (1986), 
we will label quantities that are related to the circumsolar disk with the subscript 
"1", and those that are related to the circumplanetary disk with the subscript "2". 
The reader is referred to Appendix A for a list of symbols. 

I. The creation of a protosatellite "swarm" 7 will produce a disk of prograde 
and retrograde satellitesimals extending as far as circumplanetary orbits are stable 
~ Rh/2 (Safronov et al. 1986). Initially, the capture rate is dominated by inelastic 
collisions within the planet's Hill radius of similarly-sized heliocentric planetesimals 
(see Section 2.1). As the satellitesimal disk becomes more massive its volume density 
p 2 overtakes the heliocentric planetesimal density p\. Once this happens, the dom- 
inant mass capture mechanism becomes collisions between planetesimals and larger 

satellitesimals (see Section 2.3). 

7 Although this terminology was used by Safronov to describe the circumplanetary disk, we shall 

adopt the less ambiguous choice of calling swarm objects satellitesimals. 
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2. The net specific angular momentum i z of this satellitesimal swarm is likely 
to be small; collisions between prograde and retrograde satellitesimals will result in 
a close-in, prograde satellite disk. Close to the giant planet, hypervelocity collisions 
with incoming planetesimals may shatter embryo-sized or smaller objects (~ 10 2 — 10 3 
km), and alter satellite compositions, which may ultimately result in a diversity of 
outcomes (i.e., Jovian-like vs. Saturnian-like satellite systems). 

3. Collisions between outer disk satellitesimals can lead to fragmentation, accre- 
tion or removal from the outer disk to the inner disk, where this material will be 
accreted by satellite embryos. Removal of material from the outer disk is balanced by 
planetesimal fragment collisional capture (by larger circumplanetary satellitesimals). 
Here we rely on fragmentation to decrease the size of heliocentric planetesimals, which 
are easier to capture. 

4. In the outer disk we expect that numerous ~ 1 km objects 8 (of which irregular 

satellites may be surviving members) contribute a low optical depth t z2 ~ 1CT 5 , which 

yields a surface density a 2 ~ 1 g cm' 2 and volume density p 2 ~ 1.4 x lCr 14 (<7i/crMM) 

g cm" 3 , where <j\ = o"mm = 3.3 g cm" 2 corresponds to the solids surface density of 

a minimum mass solar nebula (MMSN; Hayashi 1981). This means that even for the 
8 However, later we obtain solutions such that the satellitesimal population is made up of rubble. 

In general, satellitesimals in the meter to kilometer size range are used to capture planetesimal 

collisional fragments in the same size range. 
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case of a solids MMSN, <j\ > a 2 but typically p\ < p 2 , as noted in 1. To avoid fast 
satellite formation the mass in the extended disk of solids at any one time (excluding 
satellite embryos) must be a small fraction of a Galilean satellite (~ 10 26 g). This can 
lead to a long 10 5 — 10 6 years satellite formation time (presumably consistent with 
Callisto's partially differentiated state) controlled by planetesimal feeding. 

5. As the giant planet clears its feeding zone in a ~ 10 5 year timescale small 
planetesimals 9 are fed from the outer regions of the feeding zone and deliver enough 
angular momentum to the disk to form the regular satellites at their observed loca- 
tions. 

Satellitesimal disk collapse time: Given typical satellitesimals of size ~ 1 km, and 

a total disk mass a fraction ~ 0.1 of a Galilean satellite (which corresponds to a 2 ~ 1 

g cm" 2 ), the collapse time of the disk is of order r coH ~ P 2 /t z2 ~ 10 5 years, where 

P 2 is a characteristic orbital period of satellitesimals, and t z2 ~ 10~ 5 is the optical 

depth contributed by the larger particles (see Section 2.3). The outer disk provides 

materials to the inner disk out of which the regular satellites form. Disk material 

must then be replenished to the outer disk by planetesimal capture, and one might 

need N c = T acc jr co u ~ 10 such collisional cycles to complete the accretion of the 
9 This population of objects deserve special treatment since they will have smaller inclinations 

than eccentricities. Ohtsuki and Ida (1998) discuss the effects of small planetesimals on planetary 

spin rates. 
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Galilean satellites in r acc ~ 10 5 — 10 6 years, which may be long enough to explain 
Callisto's partially differentiated state. 

Planetesimal capture: In order to show that a GPPC model is viable one needs to 
address whether such a disk would be able to capture a sufficient amount of mass to 
form the Galilean satellites. Safronov et al. (1986) considered the following delivery 
mechanisms: 

• I. Collisions between planetesimals within R H , which may have led to the for- 
mation of a circumplanetary disk. 

• II. Collisions between large satellitesimals and smaller planetesimals, which may 
result in capture of the planetesimal and a mass inflow rate, 1^, and collisions 
between planetesimals and many smaller satellitesimals, which may result in 
capture of the planetesimal (if it encounters a mass comparable to itself) and a 
mass inflow rate I sp . 

• III. Erosion or delivery of ejecta from a collision of a planetesimal with a smaller 
satellitesimals, which may provide a mass inflow rate, I e . 

The first method of mass delivery (I) is addressed in the formation of the circumplan- 
etary disk in Section 2.1. As indicated, this mechanism is dominant until the solids 
volume density in the satellitesimal disk p 2 exceeds that of the solar nebula, p±. In 
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our model mechanism (II) determines whether enough mass can be delivered to the 
system in the planetesimal feeding time. Although mechanism (III) may provide a 
source of material, the amount of mass captured by collisional erosion is uncertain 10 , 
so we do not use I e or address it any further. We can estimate the mass delivered to 
the disk on a timescale of r acc ~ 10 6 years using (see Section 2.3) 

i-Rd 

M = 27cr acc I(R)p(R)RdR, (1) 
JRo 

where I(R) is the total mass inflow collision rate per unit area, p(R) is a capture 
probability, and i?o and Rd define the inner and outer edges of the circumplanetary 
disk, respectively. Unless otherwise noted, we generally take R to be the radius of 
the planet Rp, and Rd = Rh/2- The inflow rate I(R) can be expressed in terms of 
the disk optical depth r z2 , solar nebula solids density p±, impact velocity of incoming 
planetesimals v±, and the planetesimal and satellitesimal cutoff radii ri jmax , r 2imax 
(Safronov et al. 1986), where (see Section 2.3) 

I(R) = I sp + lib = PiV\t z2 f 2 ' max J [t z2 - In t z2 \ , (2) 

and a mass distribution with power law exponent of q± — q 2 — 11/6 has been assumed 

for both the planetesimal and satellitesimal populations. Here, r z2 = 3cr 2 /4p s r2 imax = 
10 Note that capture of material due to erosion differs from larger planetesimals captured into the 

circumplanetary disk by smaller satellitesimals, which we treat in Section 2.3. 



17 




,max 



,max 



is the vertical optical depth 



of the satellitesimal disk (p s is the 



satellitesimal density), which is assumed to be set by the size of the largest satel- 
litesimal r2 t max (Safronov et al. 1986). Ruskol (1975) gives a functional form for 
the probability distribution p(R) which measures the probability that a collision will 
result in capture (as supposed to fragmenting and scattering into escape orbits; see 
Section 2.1) 



where 9\ is the planetesimals Safronov parameter, and / = R/Rp. We can then 
evaluate the mass delivered by planetesimals primarily striking larger or comparable- 
sized satellitesimals and obtain 



This is comparable to a Galilean satellite mass for reasonable choices of parameters 
(but see Section 2.3 for further discussion). 

However, an outstanding issue is whether planetesimals may be fed into the system 

over a timescale > 10 5 years. This is because the giant planet may clear planetesimals 
11 This definition of the optical depth is only true if all the particles in the disk are of size r 2ym ax- 

For this reason, r Z 2 should be considered more a probability of collision with a satellitesimal of size 

r2,max than an optical depth. 




(3) 
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from its feeding zone on a shorter timescale. The formation of this gap was not taken 
into account in the above estimate of the planetesimal capture rate. This can be 
seen in the following simple way. A reasonable estimate for the total amount of 
mass in planetesimals in the feeding zone of the giant planets is ~ 10M e ~ 10 29 g. 
Suppose we assume that every planetesimal and satellitesimal has a size of 1 km, 
and that inelastic collisions between them always lead to planetesimal capture into 
circumplanetary orbit. Then we can ask how many times n p will a planetesimal cross 
the circumplanetary disk within -R///2 before being scattered by the giant planet. 
Ignoring gravitational focusing and taking the feeding zone of the giant planet to be 
of width A\^3Rh, we obtain 



Given an optical depth of r z2 ~ 10 5 the total amount of material captured is ~ 10 
g. But this material would be fed in a timescale 



which may be too short given our chosen system constraints. Also, since it is unre- 
alistic to assume perfect capture (i.e., p(R) < 1) and accretion efficiencies, one may 
have to replenish the feeding zone of the giant planet with planetesimal fragments 
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brought in by gas drag and enhance the solid content of the nebula with respect to 
the minimum mass solar nebula. Although we discuss this issue in Section 3, further 
work will need to be done to assess whether this is feasible. 

Angular momentum of the satellitesimal disk: A key issue is whether the deliv- 
ery of planetesimals provides sufficient angular momentum to the disk to account 
for the angular momentum budget of the regular satellites of the giant planets. As- 
suming that a significant fraction of the mass is in small planetesimals r± < 1 km 
with large Safronov parameter 9± » 1 (so that p(R) < 1), we can calculate the 
amount of angular momentum delivered to the satellitesimal disk as the giant planet 
clears its feeding zone. A rough estimate of the angular momentum budget may be 
obtained assuming a Rayleigh distribution of planetesimal velocities and neglecting 
the gravitational effects of the planet. The specific angular momentum delivered by 
planetesimals is given by (see Section 2.4) 



and the mass inflow rate is 



M = 2n X(R)p(R)RdR (8) 
where Xc(TZ) is the angular momentum delivery rate per unit area given perfect 




(7) 
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angular momentum deposition efficiency e = 1, and T(R) is the mass collision inflow 
rate per unit area. This gives a net specific angular momentum that is too small to 
create a disk of sufficient size to make the satellites unless most of the planetesimals 
originate from the outer regions of the planet's feeding zone. If we use a distribution 
of semi-major axes depleted inside an annulus of half-width ~ IRh centered on the 
giant planet, and a velocity dispersion of planetesimals v± ~ 6 km s _1 (#i ~ 100), we 
find (see Section 2.4) 



where R co u is the centrifugal collapse radius, and Mj is the mass of Jupiter. This 
suggests that it may be possible for the extended disk of prograde and retrograde 
satellitesimals to collapse to the observed Galilean satellite system size. 



Turbulence may act both to remove gas from the subnebula and to replenish it by 
feeding gas through the giant planet's gap. However, we do not rely on this gas inflow 
as a source of material to make satellites for several reasons. First, the inflow rate is 
nearly unconstrained beyond the fact that it must be very weak if ice-rich satellites are 
to survive (orders of magnitude slower than in Lubow et al. 1999; see, e.g., Makalkin 





2 Planetesimal Model 
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et al. 1999, Canup and Ward 2002). Second, the high-specific angular momentum of 
this component would lead to the formation of an extended disk around the planet 12 , 
whereas regular satellites are observed to be close-in. Third, the solid content of this 
gas component is unknown, and a significant fraction of this material might accrete 
onto the planet along with the gas. Fourth, most of the mass of solids in the planetary 
disk is likely to be in objects larger than the gas coupling size 13 , which means that 
planetesimal dynamics should be treated. 

2.1 Formation of the Circumplanetary Satellitesimal Swarm 

Planetesimal circumplanetary collisional (or collisionless) capture can be investigated 

as a three-body problem with planetesimals treated as massless particles. Inelastic 

collisions of heliocentric planetesimals within the Hill sphere of the giant planet can 

lead to capture of solids because inelastic collisions decrease relative velocities. In the 

3-body problem one can define the origin to lie at the center of mass of the system, 

and write down the Jacobi constant as 

12 It is important to note that slower inflow rates correspond to larger gaps and higher specific 

angular momentum for the inflow. Moreover, turbulence would cause the disk to spread out. 

13 Gas drag may protect sufficiently small particles, which may grow instead of fragmenting. It 

could be argued that the collisional cascade will proceed to even smaller particles in the presence of 

turbulence, but one would then need to explain how such objects formed in the first place. 
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Cj = nj(x 2 + y 2 ) + 



2GM, 



© 



2GM P 
R 




(10) 



+ 



a 



Inelastic collisions can increase this quantity (by decreasing v 2 ) to values indicative 
of a circumplanetary orbit, which allows for the possibility of planetesimal capture by 
the planet's gravitational field. Here the x-axis is directed along the line connecting 
the finite masses, a is the planetesimals distance from the Sun, R its distance from 
the giant planet, and fl\ is the uniform angular velocity of the system. Planetesimal 
collisions are inelastic because there is loss of kinetic energy through a number of 
mechanisms such as fragmentation, heating, melting, evaporation, and the like. Thus, 
collisions can occur that lead to a change in velocity < \Av\ < v such that the 
variation in the Jacobi constant (constant absent a collision) is ACj ~ — 2vAv > 0. 

We can obtain an analytical estimate of this change in energy by assuming that 
similar-sized particles collide while travelling at similar speeds v ~ v w — 
Scaling the system to units such that the gravitational constant G — 1, distances are 
in AU, masses are in terms of the mass ratio of secondary to primary /i = Mp/M Q , 
and time in terms of the orbital period (1 AU) over 2ir, then in the case of Jupiter 
one finds 



ACj w v 2 w 2 x 10 : 




(11) 
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We can compare this to the difference in the Jacobi constant of circular orbits at 
30i?j (a representative size of the Jovian satellite system) with respect to the LI 
Lagrange point ACj ~ 0.04. Thus, inelastic collisions likely allow for the capture of 
planetesimals for Jacobi constants close to the value at LI into orbits like the Galilean 
satellites with capture into orbits at further distances even more energetically feasible. 

However, mass differences between the colliding particles, as well as their initial 
energies leads us to introduce a capture probability. Ruskol (1975) derives a heuristic 
capture probability that applies to collisions between similar-sized particles, which 
she defines as those particles whose masses do not differ by more than a factor of two. 
Incoming planetesimals are assumed to have speed v w at infinity and as they approach 



the giant planet they reach a speed Vi = a/^io + v ti where v e = ^JlGMpjR. The 
velocity vectors at infinity and at distance R are randomly oriented. For a collision 
taking place at R, loss of energy from the collision allows for capture of the colliding 
bodies or their collisional fragments into circumplanetary orbits if they retain a kinetic 
energy fraction f c < (f e /fi) 2 , or more simply if their speed after the collision is less 
than the escape speed at distance R. The condition for capture is (Ruskol 1975) 

This indicates that the allowable fraction of post-collisional kinetic energy decreases 
with distance from the planet. Assuming that the probability of capture is propor- 
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tional to this term, we can write 

p(R)=Pf c (l-R/RH), (13) 

where the parameter f3 < 1 accounts for the average fractional amount of material 
captured in the collision. For our model, we assume (5 = 1/2; however, this parameter 
remains uncertain and depends on physical processes which we do not address here. 
According to this heuristic expression, the probability of capture is highest close to 
the planet and vanishes at the planet's Hill radius 14 . 

We can now estimate the amount of mass that can be delivered to the circumplan- 
etary disk via "free-free" collisions (i.e., collisions between heliocentric planetesimals 
occuring within the Hill radius of the planet, as supposed to "free-bound" collisions 
between heliocentric planetesimals and circumplanetary satellitesimals). As pointed 
out by Ruskol (1975), inelastic collisions between heliocentric planetesimals of dis- 
parate masses will not contribute significantly to the delivery of mass to the circum- 
planetary disk. Hence, one need only consider collisions between planetesimals of 

comparable mass in the range m/x < m < x m i where x ~ 2. For the planetesimal 
14 We note that it may be energetically possible for collisions between comparable-sized planetes- 
imals just outside the Hill sphere to be captured. However, in this paper, we restrict ourselves to 
the condition that only collisions within the Hill sphere can lead to capture as explicitcly assumed 
in the heuristic expression in Eq. (13). 
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mass spectrum, we use a power law of the form n{m) = Cm q , where the constant 
C can be determined by the condition for the volume density of planetesimals: 



P 



I 



n(m)mdm. 



(14) 



For the case of the circumsolar disk we have 



Ci = pi x < 



2-gi 

m l,mL- m i, 



for gi ^2; 



(15) 



[In (m limax /m limin )] for q l = 2. 



where mi : „ n and mi imaa; are the lower and upper cutoffs of the planetesimal mass 
distribution. An analogous expression can be derived for the satellitesimal distri- 
bution. One might expect that at the time of giant planet formation, most of the 
mass of solids will be in planetesimals in the 1 — 10 km size range (Weidenschilling 
1997). This may change significantly as the solids in the solar nebula evolve due to 
collisional grinding following giant planet formation, but we still use this size range 
as a reference (see Section 2.3 for more discussion). At any rate, it is fair to say that 
smaller planetesimals make our model more probable. 

As the lower limit, we take m^ min to be the gas decoupling size at the location of 
the giant planet. This is the size for which £L\t s ~ 1, where t s is the stopping time. 
The drag force on a particle depends on its size r and the drag regime. The mean free 
path of a gas molecule is \q ~ (a/1 AU) 11 / 4 cm for the standard minimum mass solar 
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nebula model (e.g., Hayashi 1981), which gives Ac ~ 1 m at the location of Jupiter. 
The stopping time is t s = p s r/pc (Weidenschilling 1977), where p and c are the local 
gas density and sound speed. At the location of Jupiter, the decoupling size is ~ 1 
m in the Epstein flow regime. 

The probability of a planetesimal collision per unit length is given by 



y(m) 



/ ' 



^nr' ni(m') dm' 



Aixp & 



2/3 ^Pi(2-gi) m 5 ^ 

— F7^— i^r, (16) 

qi - 5/d m x ^ ax - m^Zn 



which is valid for q 1 > 5/3. In this expression, the cross-section factor £ 2 for 
comparable-sized bodies, and 71 = ^1-5/3 _ ^5/3-?i_ Given Eq. (15), the rate of 
planetesimal mass undergoing collisions per unit volume can be expressed as (Ruskol 
1975) 



Ji 



mi,max 3 „2 V ^ 

y(m)vini(m)mdm = — — 7^ 

4p s ri,ef f 



The effective radius ?"i, e // is given by 



(17) 



r l,e// = < 



(<?l-5/3)(ll-6<7l) ( r l,mlx V l min) 2 
l.max l,mm 

n f r V2 1/2 s 2 

£ In (^Lmas/^brnm) 

(<?1— 5/3) (6^1 — 11) ( r l,max^l,min) 6g 11 ( r i, ma x — r l.min ) 2 



3?(2-gi)2 



ti 9l -n_ b 9l -ii 

1 .max 1 .rmn 



qi < 11/6; 

gi = 11/6; 
Ql > 11/6. 



(18) 
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Depending on the value of the power law exponent q±, ri >e // can be determined by 
the maximum particle size (q\ < 11/6), or the minimum particle size (q\ > 11/6). In 
the case qi = 11/6, most of the mass (a fraction f m — 1 — (r/r ljTnax ) 3< - 2 ~ gi ^ ~ 70%) 
lies in the top decade of the mass distribution (e.g., ~ 1 — 10 km for r^ max = 10 km). 

The flux of material T\ captured into the circumplanetary disk per unit volume 
per unit time due to free-free collisions is then the product of Eq. (13) and Eq. (17). 
For qi = 11/6, the mass delivery rate per unit volume can then be written as 

Tx = Pi(R)Ji = PiPi^i7i77- hi ( Tl ' max j , (19) 

which has been written in terms of the local vertical optical depth of solids in the 
nebula r z \ = 3<xi/4p s ri (see footnote 9), where o\ = p\H\ is the surface density of 
solids in the planetary zone, Hi = \/2vi/Q P\V\/A is the local scale-height of solids 
in the nebula, and P 1 is the planet's orbital period. A rough estimate of the amount 
of mass delivered per year to the circumplanetary disk is 



M 1 = / T\ dV w 247T7i^ —In I— J / (l — x/x H )dx 

J PlPsfl,max \fl,minj J l/20i 1 + x 

V3.3gcm- 2 y \r 1<max J 
where xh = Rh/2@iRp, Ps = 1 g cm -3 , and we assume 9\ = 4. The timescale for 
the giant planet to clear its feeding zone is At ~ 10 5 years (e.g., Tanaka and Ida 
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1997; Charnoz and Morbidelli 2003), so this is presumably the timescale in which the 
circumplanetary disk should be populated. Assuming our fiducial choices for the disk 
parameters, one finds that the total mass delivered then is At x M x ~ 10 23 g, which 
is considerably smaller than a Galilean satellite. This would seem to pose a problem; 
however, we show below that once the circumplanetary disk becomes sufficiently 
massive, capture of planetesimals by inelastic collisions with circumplanetary disk 
particles (free-bound collisions) becomes the dominant mass delivery mechanism. 

2.2 Planet esimal Capture Due to the Circumplanetary Disk 

The capture of material due to free-free collisions may form a circumplanetary disk 
around the giant planet. As the density of the circumplanetary disk p 2 increases, 
it may eventually become larger than p 1 . Once this happens, the primary delivery 
mechanism may become free-bound collisions (Ruskol 1975, Safronov et al. 1986). 

Along the lines of the prior section, except that satellitesimals are confined to a 
circumplanetary disk with scale- height H 2 , we first estimate the mass rate per unit 
area undergoing collisions between planetesimals and satellitesimals of comparable 
size. In analogy to Eq. (16), the probability of a collision of a planetesimal crossing 
the circumplanetary disk in the normal direction is given by 
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pi \ /"*% *n ( 'w ' ( 3 V /3 W2-g2) m5/3 ^ 2 

£{m) = / £?rr H 2 n 2 {m)dm = T7^72 



'm/ x V<±"^/ <?2-5/3 m^-m^' 

(21) 

which is valid for g 2 > 5/3. Here, n 2 = C 2 m~ q2 is the mass distribution of satellitesi- 
mals with the coefficient C 2 defined as in Eq. (15), 72 = x 92-5 / 3 — x 5//3_92 , and iJ 2 and 
(T2 are the scale-height and surface density of solids in the protosatellite disk, respec- 
tively. The total number of planetesimals impacting the circumplanetary disk per 
unit area over the mass interval m to m + dm is vini(m)dm. Integration over m gives 
the mass rate per unit area of collisions between planetesimals and satellitesimals of 
comparable size 



h — \ £(m)vin 1 (m)mdm = 72 ~ H 2 J 2 , (22) 

Jm 2 , min 4 p s r 2 , e ff 

The effective radius r 2e tt can be expressed in terms of the appropriate limits 



r 2 ,eff = n < 



(11/3-91-92X92-5/3) 1 . , 11/0. 

?(2-9l)(2-9 2 ) (ll-39l-392)_ (11-39X-392) HI T «/2 ^ 11 /<->, 

2,max 2,min 

92-5/3 1 _ , _ .ii/q., (23) 

3^(2-9i)(2-9 2 ) ln(r 2 , max /r 2 , mi „) q±-rq 2 - LL/O, 

(gi+92-ll/3)(92-5/3) (ra.^ra.min) 11 - 39 !- 3 ^ n ,„ s 11/0. 

C(2-9l)(2-9 2 ) (ll-3gi-3g 2 ) (H-391-392) VI T (/2 ^ 1 V <J, 

2, max 2,mm 

7? = iV 6 " 391 _ -G-^iVt- 6-3 * 2 _r 6_3 ' 32 'i (94^ 

/v — V' l,mac ' IminJV 2,max 2,min I ' 
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The choice of r 2 ^ min is problematic given that the gas surface density in the circum- 
planetary disk is left unspecified. We shall, unless otherwise noted, arbitrarily choose 
this size to be ~ 1 cm, though our results are not dependent on this parameter given 
our choices for q 2 — 11/6 and r 2tTnax ~ 1 km. 

We now consider the conditions under which capture by the circumplanetary disk 
begins to dominate the capture rate. We can express the capture rate of material due 
to collisions of comparable-sized bodies as 

Tcs = Fi + F 2 = Pi Ji + -K-J2 = 7 7i + 7. 72 , 25 

2 4 ft \r\ eff 2r 2 ,eff J 

where we have assumed that the typical amount of mass involved in a free-bound 
collision is roughly half that of the free-free case (Ruskol 1975). This is because we 
are only considering comparable-sized collisions for both capture rates. From Eq. 
(25), we see that in the case of collisions between particles of comparable size the 
condition that the dominant capture mechanism involves the circumplanetary disk is 

> 9 7l r 2,eff _ ( r 2,max \ ^ In (^l.maa/^l.mm) (r)a ^ 
92 £ *Pl ^Pl 7— r ~ Pi, [Zb) 

l2n, e }f \ri, m axj in.{r 2 , max /r 2 , m in) 

where the next to last term on the right side applies for g x = q 2 = 11/6. Independent 
of the uncertainty in the lower bound of the size distribution for satellitesimals, we 
see that the circumplanetary disk dominates the capture of planetesimals when the 



31 



volume density of solids in the disk becomes comparable to that of the solar nebula 
(Safronov et al. 1986). 

To determine the implied volume density of solids in the circumplanetary disk 
that satisfies this condition, we estimate the mean surface density of solids in the 
protosatellite disk to be 



4M d 4MiAt 
° 2 = n(Rl-ARl) ~ „(R* H -AR%) ~ °-° lgCm " ' (2?) 



where we have taken the mass delivered in At = 10 5 years due to both mechanisms 
to be approximately M ± At, we have made use of Eq. (20) with 6>x = 4, and we have 
ignored collisional removal. This means that a 2 « U\. Since the critical density such 
that capture by the circumplanetary disk becomes dominant is p2, C rit = ^/^^rit = 
pi, for our disk parameters the scale- height of solids in the circumplanetary disk must 
be less than 

B^u < wo ( ( «*e) Rj , (28) 

V 3.3 gem 2 J \r limax J 
and the critical density must be greater than 

P2,cnt > 1.4 x 10- 14 ( ai - ) gem" 3 . (29) 
\3.3gcm 1 J 

We stress that this is the condition for disk capture to become dominant. Once this 
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happens the subnebula surface and volume densities may increase significantly. As 
we show later, in general a 2 < <J\ to be consistent with chosen system constraints. 

2.3 Total Mass Delivered to the Circumplanetary Disk 

We now remove the requirement that colliding particles be of comparable mass, calcu- 
late the mass inflow rate due to the entire mass spectrum, and then use this quantity 
to estimate the mass delivered to the circumplanetary disk over the planetesimal 
feeding timescale. 

Here there are two cases to consider. For a given planetesimal of mass m, we 

estimate the mass inflow collision rates per unit area (henceforth inflow rates) I sp and 

lib due to interactions with satellitesimals with w! < m, and m! > m, respectively. 

These rates include the contribution due to similar-sized bodies 15 . In the regime we 

are in, we will see that the inflow rate >> I sp since the optical depth r z2 << 1 
15 In prior sections, we estimated the mass inflow rates for two mechanisms for comparable sized 

bodies. In the case of collisions between planetesimals in the Hill radius of the giant planet, collisions 

between comparable sized objects are considered because these incur the largest velocity changes 

resulting in an increased likelihood of capture. In the case of collisions between planetesimals and 

comparable-sized satellitesimals, the mass inflow rate for q\ = q 2 = 11/6 using Eq. (22)-(24) is 

t nr. P\V\(J<2, . ( T2^max\ „„_ / f2.max\ ^ , / ^2, max \ 

Ics ^ 0.2— — H In I — : ~ 0.25piwir z2 — In — , 

Ps[J , l,max'^2,max ) \ f*2,min J \^l,max J \ ^2, mm / 

which is in agreement with Eq. (6) in Safronov et al. (1986). 
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given the constraints of the problem. 

For an incoming planetesimal of mass m to stand a chance p(R) of being captured 
by satellitesimals smaller than m, the condition is that nr 2 H c p c > m (Safronov 
et al. 1986), which simply says that in order for the mass m to have a chance 
of being captured into the circumplanetary disk by particles smaller than itself, it 
must encounter a mass comparable or larger than itself when crossing the disk. The 
characteristic size for which all particles with r < r c stand a chance of being captured 
into the satellitesimal disk is 

r c /r Xmax = [r z2 { 1 - (r Xmm /r c f-^ }] , (30) 

where we have taken (cf. Eq. [15]) p c = C^m 2 " 92 — rnl~J^ n )/(2 — q 2 ). If we assume 
for simplicity that r c » r 2im in, then for our fiducial case of q 2 = 11/6 one has 
r c ~ Tz 2 r2,max, which suggests that unless the disk solids surface density <t 2 is very 
large, only the smallest particles are assured of capture. For a given q 2 in which 
most of the mass is in the largest particles, the inflow collision rate is small (i.e., for 
Tz2 << 1)- Under these circumstances, there is a probability that all bodies below 
this size are captured into the disk (so that £(m) = 1). The inflow rate on the 
circumplanetary disk is then just given by an integral over the number of impacting 
bodies on the disk's surface 
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I sp = / vini(m)mdm ~ P\V\T^ h I ' maa: I = p\V X r z2 I '"^ I , (31) 

JO \^'l,max J \ ^l,max / 

where q± = q 2 = 11/6, and we have made use of Eq. (30) in integrating from 
to m c . For a sufficiently broad size distribution for the satellitesimals, 
(see footnote [12]). Given that in order to satisfy the long timescale of formation of 
Callisto r z2 ~ 10~ 5 , I sp is likely to be small as indicated above. We can see this if we 
examine a case such that r z2 = 1, which means r c = r 2) max- Then all planetesimals 
with r < r 2j max could be captured into the disk. However, this condition is quite 
restrictive. For instance, a disk with the mass of the Galilean satellites out to ~ 30Rj 
has <j 2 ~ 10 3 g cm -2 , which implies r 2 ^ max ~ 1 - 10 m. Even if such a disk could 
be produced, the formation time for the satellites would be short, which would be 
inconsistent with Callisto's partially differentiated state. 

Now we consider the mass delivery to the circumplanetary disk due to the inter- 
action of planetesimal fragments with larger satellitesimals. The expected number 
of collisions of a disk crossing planetesimal of mass m with satellitesimals of mass 
m' > m is 



S(m; m > m) 
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which is Eq. (11) of Safronov et al. (1986), where we have corrected an error in 
one of the exponents in that expression. Thus, the frequency of interaction of bodies 
with m < m 2jma x is proportional to m 5 / 3 ~ 92 (Ruskol 1975). Integrating from m 2>min 
to m 2 , max we find the three cases of interest for the mass inflow rate of planetesimals 
colliding with satellitesimals of larger mass 
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where the lower limit of integration has been taken to be r c . For our chosen values 



of qi — q 2 — 11/6, the inflow rate due to larger bodies is 



fr 2max \ 112 

hb ^ PiV\T z2 I J [t z2 + In r" 1 - l] , (35) 

so that the total inflow rate into the circumplanetary disk for qi = q 2 = 11/6 is 



I — I S p + I lb ~ P\V\T z2 



r 2. 



Tl,max 
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[r z2 ~ In r z2 ] . 



(36) 



We see from this expression that for a low optical depth disk the mass delivery of 
incoming planetesimals by larger satellitesimals is more efficient than delivery by 
smaller or comparable sized bodies, as already mentioned. Note that in Eq. (36), the 
second term becomes zero for values of r z2 > 1 since this corresponds to the situation 
in which r c ~ r 2ima:r {i.e., all bodies smaller than r 2jm , a:E could be captured into the 
disk). 

The mass accretion rate obtained by integrating Eq. (36) over the disk area is 



where p(R) is given by Eq. (13), I = R/Rp, and r 2 and r\ are the satellitesimal and 
planetesimal size distribution cutoffs. We have expressed the quantity p^vi in terms 
of the surface density of solids in the nebula a±, and the planetary period Pi. 

We can now obtain solutions for the mass captured by the circumplanetary disk 
assuming that the mass delivery rate to the disk is balanced by the mass removal 
rate to the inner disk (and used to form close-in, prograde regular satellites). The 
mass removal rate is given by M ~ M sw /t coUi where M sw is the steady state mass of 
the satellitesimal swarm (which excludes satellite embryos), and the collisional time 
is t co u ~ P 2 /r 22 . We choose the characteristic orbital period in the disk to be P 2 ~ 1 
year, which is the orbital period at R ~ Ph/4 (where the maximum of p(R) x R 



M 




(37) 



37 



is located for 9\ ~ 75). Given that the total mass to be delivered in the satellite's 
accretion time r acc is M ~ 10 26 g, we must have M = M sw + M x r acc . This simply 
means that the disk mass must satisfy M sw = M/(l + N c ), where N c = T acc /r co u is 
the number of collisional cycles. We can also write this equation in the form 

N 

Mr acc - -^W M = °> ( 38 ) 

which allows us to express r 2 , t z2 , and a 2 in terms of N c , so that equations (37-38) 
can be solved by iteration as functions of N c , for specific choices of the solar nebula 
parameters r±, 9i and U\. Here we assume that the optical depth and solids surface 
density of the circumplanetary disk are constant not only in time but with position 
as well. 

Figures 1 and 2 correspond to a case such that the satellite accretion time (which 
is set by the planetesimal feeding time) is r acc = 10 6 years, the solids surface density 
in the solar nebula is ten times MMSN, and the Safronov parameter 9\ — 4 (v w = 21 
km s _1 ). In Fig. 1 we plot the satellitesimal cutoff size r 2 (solid line) and surface Figure 1 
density cr 2 (dotted line) as a function of r\ that satisfy the system constraints. In 
this case solutions for which N c ~ 1 — 100 correspond to the range of r\ shown. 
Larger values of r 2 correspond to larger t coU (smaller t z2 ) and hence fewer collisional 
cycles, as indicated in Fig. 1 and by the dotted curve in Fig. 2. The solids Figure 2 
surface density a 2 changes only by an order of magnitude over the same range in 
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which the satellitesimal maximum cutoff size varies by a factor of ~ 1000. The 
greater variation of r 2 compared to a 2 is a result of the of the stronger dependence 
of r 2 oc A^~ 1 (l + iVc) -1 compared to a 2 oc (1 + iV c ) -1 . It should be pointed out 
that there are solutions for N c ~ 1. This corresponds to a case in which the steady 
state disk mass is ~ 0.5 the mass the Galilean satellites. It is unclear, however, 
whether the required size of r 2 ~ 10 km (for 6\ = 4) is too large to be realistic. 
We do not attempt to model their size distribution, but a population of such large 
satellitesimals may be difficult to maintain given fragmentation from hypervelocity 
collisions from Roche-lobe interlopers. We discuss this point in more detail at the end 
of this section. Here we note that the range of possible values of r 2 is quite broad, 
so that an independent constraint for this quantity is essential. It is interesting to 
note that there are solutions such that the circumplanetary population is composed of 
rubble r 2 < 10 m. This can work because for these cases the optical depth of the disk 
is relatively high, so that enough mass may be captured from the small mass fraction 
at the low end of the incoming planetesimal mass spectrum (given that r\ » r 2 ). 

Next we do a case in which <7i is two times the solids MMSN, r acc = 10 6 years, 
and we assume Q\ = 75 (v w = 5 km s _1 ), implying a fairly "cold" population of 
heliocentric objects (Figs. 3 and 4). We find solutions for ri ~ 5 — 10 km, which Figure 3 
correspond to N c in the range 5 — 100. In particular, a solution exists where Figure 4 
the satellitesimal cutoff size is r\ — 1 km for a planetesimal cutoff size of r\ = 10 
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km. In this case, t Z 2 ~ 5 x 1CT 6 implying that timescale for the collisional cycle is 
Tcou ~ 2 x 10 5 years, so that N c ~ 5 collisional cycles are required to complete the 
formation of the satellites. We note that the range of variation of both r 2 and <r 2 is 
significantly smaller for this case, which suggests that large d\ cases may be easier to 
constrain. 

To understand better the effect of the velocity dispersion of the heliocentric plan- 
etesimal population on the mass delivery rate, we plot in Fig. 5 the total mass Figure 5 
delivered to the circumplanetary disk over and accretional time of r acc = 10 6 years as 
a function of the Safronov parameter 9 1 . In this case we fix r\ and r 2 at 10 and 1 km 
respectively. All of the cases presented assume that t z2 = 5 x 10~ 6 so that N c = 5. 
The solid curve corresponds to a solids MMSN (% M ), the long dashed line corre- 
sponds to a case of <j\ = 2<7mm, the dotted line to <xi = IOomm, and the short-dashed 
line to <Ti = O.Iomm- From this plot we can see that for a choice of 9\ = 4 to have 
a solution for r x = 10 km one would require a case with greater nebula solid content 
than ten times MMSN (cf. Fig. 1). Conversely, one could have a slightly colder 
population (9i ~ 8) that satisfy the constraints for the case of ten times MMSN. 

Finally, in Fig. 6 we plot the mass captured by the circumplanetary disk in one Figure 6 
collision time for r acc = 10 6 years as a function of collisional cycles iV c . This mass 
simply corresponds to the steady-state disk mass M sw as it should. All prior solutions 
for a range of r\ lie on this curve. For instance, the case for twice MMSN and r\ = 10 
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km for which N c ~ 5 collisional cycles implies that M ~ 0.2 Galilean masses are 
captured per collision time t coU ~ 2 x 10 5 years (corresponding to M ~ 5 x 10 20 
g yr^ 1 ). This means that at any given time the total satellitesimal disk mass M sw 
(excluding embryos) is ~ 20% of the total mass of the satellites. 

An important caveat is that all these cases use accretion times that are signifi- 
cantly longer that feeding zone clearing timescale, which means that we are assuming 
some mechanism to replenish it. It is also important to point out that at these late 
times modeling of the planetesimal mass spectrum would involve different regimes 
with different power law exponents (see, e.g., Charnoz and Morbidelli 2003, Kenyon 
and Bromley 2004, Pan and Sari 2004). In the next section we consider small values 
for 7~i ~ 1 km to indicate that following giant planet formation a significant fraction of 
the planetesimal mass in the feeding zone of the giant planet may cascade back down 
to small sizes. Although our solutions so far have involved somewhat larger values of 
ri, in general only planetesimal fragments < 1 km may be captured and can be used 
for satellite formation. One might expect even more mass in smaller particles, which 
as we point out below would benefit our model. 

Finally, one would still need to solve self-consistently for the size of r 2 . This size 
might be set by a competition between the growth time of satellitesimals versus the 
time it takes for catastrophic collisions with incoming planetesimals to disrupt them. 
The growth time of r 2 is ~ t co u which may be, as for the cases in Figs. 3 and 4, 
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quite long ~ 10 5 years. We can estimate the lifetime of a satellitesimal of size r 2 as 
the time it takes for a bigger planetesimal to collide with it. Given a planetesimal 
size distribution with upper cutoff of r\ = 10 km and power law exponent q\ = 11/6, 
the time for a large planetesimal to shatter a ~ 1 km satellitesimal can be estimated 
to be T frag ~ 10- 3 Pip s (rir 2 ) 1/2 /ai ~ 10 3 years for a solids MMSN. This is ~ 100 
times faster than its growth time, suggesting that the cutoff in r 2 should be smaller 
than 1 km, and invalidating the assumptions used to obtain model solutions. This 
result is a direct consequence of the model constraints coupled with a small capture 
probability p(R) « 1. It may be possible to circumvent this problem by speeding 
up the accretion of satellitesimals using circumplanetary debris, and increasing the 
capture probability by using a large value of the Safronov parameter 6\. It is also 
possible that a different process, such as the capture of planetesimal fragments or 
a gravitational instability, is responsible for setting the satellitesimal cut-off size, or 
that that the satellitesimals are composed of collisional rubble, as indicated by some 
of our solutions. However, we leave investigation of these possibilities for further 
work. 

So far we have ignored the issue of the angular momentum of the satellitesimal 
swarm, and found solutions subject to the constraint that enough mass be delivered 
to form the regular satellites absent a gap in the planetesimal disk. We obtained 
solutions such that r\ > r 2 ~ 1 km as well as solutions where r\ » r 2 ~ 1 m. 
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This is because lower values of T\ would deliver too much mass. However, given the 
uncertainty in the calculation this can not be taken to be a binding conclusion. Below 
we attempt to match both the mass and angular momentum of the satellites. In this 
case, we need to deliver enough mass to make the satellites as the giant planet opens 
a gap in the planetesimal distribution. For this reason, next we consider cases such 
that r 2 > r\ and Q\ » 1, both of which facilitate capture. Allowing for larger values 
of r\ would require us to use more planetesimal disk mass by an amount that depends 
on the (unknown) value of q. For simplicity, below we neglect these complications. 

2.4 Angular Momentum Delivery 

In this section we show that it may be possible to deliver enough angular momentum 
to the circumplanetary disk to explain the observed angular momentum budget of 
the regular satellite systems. At first this might appear to be a problem given that 
the circumplanetary disk is composed of both prograde and retrograde satellitesimals 
with small net angular momentum. However, as the giant planet clears planetesimals 
from its feeding zone on a timescale comparable to the satellite formation timescale, 
inhomogeneities in the planetesimal disk can pump up the angular momentum de- 
livered to the circumplanetary disk (a point first made by Lissauer and Kary 1991 
in the context of planetary spins), which may account for the observed radial extent 
of the satellite systems provided that most of the solids that make up the satellites 
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originate from the outer regions of the giant planet's feeding zone. 

Here we assume that the bulk of the angular momentum that goes into the regular 
satellites arrives late as the planet clears its feeding zone, and that the planetesimal 
population at this time is composed of small particles r\ < r 2 < 1 km with a velocity 
dispersion characteristic of a fairly cold population (#1 ~ 100, but below we still use 
the "high" dispersion limit) with lower inclinations than eccentricities (v z = 0.1t>i) 16 . 
This implies that the capture probability p(R) < 1. Furthermore, we assume that on 
average a free-bound collision deposits a significant fraction e ~ 0.5 of the angular 
momentum of the interloper. This is because colliding planetesimals may deposit 
angular momentum in the circumplanetary disk even if they (or their fragments) are 
not captured 17 . 

We provide an estimate for the amount of angular momentum delivered to the 
circumplanetary disk using a method which is related to planetary spin calculations 
(e.g., Lissauer and Kary 1991; Dones and Tremaine 1993, hereafter DT93). We 
calculate the rate of accretion of mass and angular momentum by a (2-D) disk of 
satellitesimals. Our simplifying assumptions are that the planet is on a circular orbit, 
16 Strictly speaking this implies that v\ < Q\Rh- 

17 One might consider the case e = p(R) as a lower bound on the angular momentum deposition 
efficiency. This would require that a planetesimal (or its fragments) be captured in order to deliver 
its angular momentum to the disk. For our choice of parameters p(R) < e. 
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that incoming planetesimals do not interact with each other, and that o\ is uniform 
over the planet's (undepleted portion of the) feeding zone. Incoming particle velocities 
increase as they approach the planet. This prevents us from considering significantly 
larger values than 9 1 ~ 100, and leads us to adjust the velocity dispersion using 



vi = a/v 10 + , with v 10 and v e as previously defined. 18 In addition, the giant planet's 
mass and the properties of the circumplanetary disk are assumed to be constant. With 
these assumptions, the problem can be cast in terms of two dimensionless parameters 

where b is the ratio of radial distance in the circumplanetary disk to the Hill radius, 
and s measures the importance of the differential rotation of the circumsolar disk 
relative to the velocities of incoming planetesimals (DT93). The parameter b is a 
measure of the importance of the planet's gravity relative to that of the Sun at R in 
the satellitesimal disk. For our problem, 9\ > 100 implies v\ < 6 km s~ x and thus 
s < 6. At Rh/4: (where p(R) x R peaks for Q\ ~ 100), b/s < 1 and bs 2 > 1, implying 
we are in a regime in which dispersion is still high and gravity is somewhat important. 

The accretion of mass and angular momentum by the circumplanetary disk can 
18 It should be noted that the satellitesimals Keplerian velocity is not much smaller than v\. This 

could affect the calculation for the net angular momentum deposition. We leave investigation of this 

complication for further work. 
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be described by a 3-D distribution $(R, £ z ), where $(-R, £ z )d£ z is the collisional inflow 
rate per unit area of mass with specific angular momentum £ z in the range £ z and 
£ z + d£ z . For our problem, we have 



where the azimuthal angle is measured from the origin of the system centered on the 
planet such that (x, y) = i?(cos 0, sin 0), the x — y plane is the planetary orbit plane 
(which coincides with the circumplanetary disk plane), the x-axis points radially 
outward, the y-axis points in the direction of orbital motion, 8 is the Dirac delta 
function, and H(—z) is the Heaviside step function, which selects only planetesimals 
with velocities directed towards the disk plane. The notation z — > + is meant to 
indicate that we are approaching the disk plane from above (z > 0). The probability 
of a free-bound collision in crossing the disk t($) is given by 



where v r — x cos + y sin is the radial velocity, is the azimuthal velocity, and we 
have assumed that the planetesimal radial velocities dominate their motion (which 
implies a shallow angle of entry). In Eq. (40), f is a Schwarzchild distribution 




(40) 



f (x, y,z -> + , x, y, z)5(£ z - xy + yx - fi!i? 2 )|i|if(-i), 




(41) 
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function, describing the mass of a uniform disk of particles in a three-dimensional 
phase-space volume element dxdydzdxdydz. 



fo(x,y,z,x,y,z) = 

(42) 



Po 
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Qjz 2 x 2 (y + T^x) 2 z 2 



2v 2 2v\ 1v\{\ -\Y) 2v 2 _ 
which may apply sufficiently far from the planet. In this expression, p — \P^Ho\jP\v z 
is the midplane planetesimal volume density, and we have assumed that the x- 
component of the velocity dispersion v x = v\, and v y = vi(l — ^T) 1 / 2 , where T = 
— (dlnfl/dlnR) = 3/2 for a Keplerian disk. 

Integrating Eq. (40) over z, defining the velocities (x,y) = (vcosip, vsinip), and 
replacing the azimuthal angle <fi with x = 7r + — ^ so that v r = —v cos x, we can 
write in polar coordinates (DT93) 



$(R,£ Z ) = / #/ cos X dx 

[2^{l-\Y)\ 1,2 vlv z pJo J-,/2 

i 

poo 

/ v 2 e-^5(£ z - Rvsm X -ttiR 2 )dv 

J\£ z -niR' 2 \/R 



(43) 



and 



_ v 2 cos 2 ip [v sinip — TfliRcos (x + ip)] 2 

+ via - §f) • m 
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Generally, if the velocity dispersion is large enough or for small b, the quantity 
b/s « 1, and Eq. (43) can be simplified by expanding Eq. (44) to first order in b/s. 
However, for our problem b/s > 0.1 so that we expand Eq. (44) to fourth order to 
obtain a better estimate. Carrying out the expansion and integrating over \ an d v 
(see Appendix B) 

*(R, «) = r ; r -TJ-J K ^) + K x {g)0L + K 2 (g)a 2 + 

R 2tt(1 - |r) ] ' v z Pi Jo 9W 



R [2tt 

K,{g)c? + K 4 (g)a 4 )e-^ + (K 5 (g)a 2 + K 6 (g)a 3 + (45) 

where the Ki are functions of the angle ip, T, and b/s. Here a = (£ z — ^liR 2 ) / y/2v 1 R, 
Ei is the exponential integral of order 1, and 



^) = cos 2 ^ + ^^. (46) 
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To obtain the mass capture rate of planetesimals with specific angular momentum 
\£ z \ > £ we define the mass inflow collision rate per unit area T{R,£ Q ) 



i poo i p-to 

1(R, e Q ) = — / $(i?, 4) d£ z + — $(R, 4) d4 = 

" (47) 
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where $(-R, a') is identical to Eq. (45), except with the replacements K\ — > —Ki, 
K 3 — > —if 3 , if 6 — > — iQ, an d — (4 + ^i-R 2 )/\/2vi-R- Integrating over the disk 
area, we get the mass capture rate 



M(£ ) = 2ir l(R,£ )p(R)RdR = 

Jro 



,. 7 , (48) 

where p(i?) is given by Eq. (13), and the M.i are given in Appendix B. The total 
mass delivered to the circumplanetary disk is then given by M = M sw + Mr acc , 
where as before M sw is the steady-state circumplanetary disk mass (excluding satellite 
embryos) . 

The angular momentum rate delivered to the circumplanetary disk L z can be 
calculated from the first moment of the distribution $(i?, £ z ). The angular momentum 
delivery rate per unit area Xl(R,£ ) assuming perfect deposition efficiency is 



i r°° i /"—to 

i L (R, ( ) = 7T £ MR, Q d£ z + — £MR, Q d£ z = 

-.2 p2 [ roc poo 

n l R 2 l(R) + ^ — <^ / a$(R, a) da- / a'<Z>(R,a') da' 

71 Ja J a' 



(49) 

which when integrated over the circumplanetary disk yields 
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f IfCT T O f 1! 

L z (£ ) = 2ir eX L (R,£ )RdR = 1 z \ * / —RdRx 

J Re (1 - |r) 1 P x J Re v * 

(50) 

where the Ci are given in Appendix B. The total angular momentum delivered to 
the circumplanetary disk is then L z = Li + L z r acc , where Li is the initial angular 
momentum of the disk, which is taken to be small. The specific angular momentum 
is (4) = L z /M. 

We integrate Eq. (48) and (50) over angle ip and radius using standard numerical 
techniques (note that, although v\ varies with R, we held the ratio v\jv z constant in 
these calculations). In Fig. 7 we plot both the mass M and z-component of angular Figure 7 
momentum L z delivered to the circumplanetary disk over the planetesimal feeding 
time r acc = 10 6 years as a function of the feeding zone gap-size R gap = a/^o/^i 
with Rq = R P and Rd = Rh/2- It is important to note that this is a practical 
definition based on the observation that for R gap > 2.5Rh the mass inflow rate drops 
to nearly zero. Because in the high dispersion weak gravity limit large magnitude 
(planetesimals may deposit positive or negative) angular momentum contributions 
are not limited to the outer regions of the feeding zone of the giant planet, this 
calculation likely overestimates the increase in specific angular momentum as a gap 
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is opened 19 . Here we are mainly interested in the relative behavior of the angular 
momentum curves to the mass curves. The mass and angular momentum are plotted 
so that unity corresponds to the total mass and angular momentum contained in 
the Galilean satellites, M sats ~ 4 x 10 26 g, L Z)Sats ~ 4 x 10 43 g cm 2 s" 1 , which 
yields (i z ) sats ~ 1 x 10 17 cm 2 s _1 . In this figure we plot curves for combinations of 
surface density o\ and angular momentum deposition efficiency e such that solutions 
matching the mass and angular momentum in the Galilean satellites were found for 
gap sizes R gap ~ IRh- This implies that for our chosen parameters the bulk of the 
solids that form the Galilean satellites can be delivered as the giant planet clears a 
gap in its feeding zone. This can work provided that a significant fraction of the mass 
of the solids is in planetesimal fragments with r\ < r 2 with low inclination orbits 
v z = O.lvi, that collisions are efficient at depositing angular momentum e ~ 0.5, and 
that the nebula is enhanced in solids with respect to MMSN. Given a value for the 
satellitesimal size cutoff of r 2 = 1 km, we then have r Z 2 ~ 10~ 5 , t co u ~ 10 5 years (since 
the dynamical time in the satellitesimal disk is ~ 1 year), N c ~ 10, M sw ~ 0.1M sats 
as before (which is the value the mass curves [dotted lines] approach as the gap size 
Rgap encompasses the giant planet's feeding zone), and Li < 0.1L sats . We can also 
define the quantity 

19 Ohtsuki and Ida (1998) argue that such a boundary should be implemented in terms of the 
Jacobi constant. 
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A = <JiT z2 (^j Tacc = a,P 2 (j^j N c , (51) 

where a± is taken to be in units of the solids MMSN (ct M m = 3.3 g cm" 2 ), and r acc 
and P 2 are in years. The pairs of curves in Fig. 7 correspond to A = 280 (e = 0.8, 
(Ti = 2.8ct M m), A = 450 (e = 0.5, a x = 4.5<t M m), and A = 1200 (e = 0.2, a x = 12ct mm ). 
If we change the values of parameters we can obtain approximately the same curves 
as before by keeping A constant. One can also consider longer feeding timescales, but 
such solutions might be unrealistic. 

In Fig. 8 we plot the specific angular momentum as a function of the gap size for Figure 8 
parameters such that the mass and angular momentum of the Galilean satellites are 
obtained for (£ z ) / (£ z ) sats = 1- As one might expect, the specific angular momentum 
curves first increase with gap size, reach a maximum, and begin to decrease as the 
gap size chokes off the mass inflow. It is interesting to note that for the case of 
nearly perfect angular momentum deposition efficiency e ~ 1 it would seem possible 
to match the constraints even without a gap. However, it remains the case that to 
satisfy the system constraints more generally most of the angular momentum must 
be delivered as the giant planet opens a gap in the planetesimal distribution. 

In Fig. 9 we plot the specific angular momentum delivered as a function of the Figure 9 
satellitesimal disk size Rd for several values of the gap size. This figure shows that 
one may need an extended disk ~ Rh/2 in order to be able to deliver sufficient mass 
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and angular momentum to form the regular satellites of Jupiter (and Saturn 20 ). 

Finally, in Figs. 10 and 11 we increase the size of the inner disk boundary Figure 10 
R and compare the relative behavior of the angular momentum and mass curves as 
we did before. Strictly speaking, this procedure would be tantamount to opening Figure 11 
a gap in the heliocentric planetesimals in the limit of small velocity dispersion and 
weak gravity, which is inconsistent with the assumptions employed to obtain these 
curves. Nevertheless, because opening a gap would preferentially decrease the inflow 
rate close to the planet, it is still instructive to consider this case. Here we simply 
re-scale the size of the gap by R gap = 5R . The results are qualitatively similar to 
before, except that the effect on (£ z ) is less pronounced and a solution could not be 
found in the case of e = 0.2. 

It is important to note that these results employ r acc = 10 6 yrs. An outstanding 

issue is whether planetesimals may be replenished by gas drag in the feeding zones of 

the giant planets on a timescale > 10 5 years. By this time the giant planet would have 

opened a more or less clean gap in the gas as well as the planetesimal disk (depending 

on the unknown turbulence parameter), which means that even meter-sized particles 

may not be well-coupled to the gas in the feeding zone of the giant planet. In this 

limit the gas density may not be responsible for selecting the sizes of planetesimals, 
20 Saturn's satellite system might require one to adapt the present formalism to the case of a giant 

planet with significant obliquity. 



53 



and the collisional cascade may proceed to possibly even smaller sizes. In regions 
outside the gas gap, resonances may capture inwardly drifting planetesimals, which 
might prevent them from replenishing the feeding zone. However, collisions or stirring 
by larger particles may free them and allow them to continue their inward migration. 
As a result, this problem is unique in that it spans both gas- free and gaseous regimes. 

Ignoring the gap in the gas, the criterion for feeding zone cleanup is (Tanaka and 
Ida 1997 and references therein) 



Ul + 0A6(T h/T K )^- 1 j / (1 + 0.46(T V^) 2/3 ) 1/4 > v/K (52) 



where T = ImjCD^f' 2 p g [hvK) is the characteristic gas drag time, m is the particle 
mass, p g is the gas density, Tk is the orbital period of the planet, h = (Mp/SMq) 1 / 3 , 
C D = 0.44, v K is the Keplerian velocity of the planet, c is the gas sound speed, and 
7] = c 2 jv\. In the above equation, it is assumed that i < e. If the criterion for 
cleanup is not satisfied, planetesimals outside the planet's orbit can enter the feeding 
zone of the planet. For Jupiter, the size that can replenish the feeding zone of the 
planet is ~ 1 m, where we have used a minimum mass solar nebula (Hayashi 1981). 
But simulations indicate that bigger sizes may be able to replenish the feeding zone of 
the giant planet (Tanaka and Ida 1997). It is also possible that planetary migration 
may help in this regard, or that one planet may stir up the planetesimal disk and 
replenish the feeding zone of its neighbor. 
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3 Discussion and Conclusions 



We have focused on outstanding issues that a planetesimal collisional capture model 
faces before it can provide a viable alternative to the model developed in MEa,b. Its 
strength is that even though it assumes that turbulence can drive the evolution of 
the subnebula disk in a shorter timescale than that of satellite formation, a turbu- 
lence parameter appears nowhere in the analysis. While several of the parameters 
employed herein are at best poorly constrained at this time, none of them are incal- 
culable, and all are based on sound physical processes. We have investigated whether 
the mass and angular momentum budget in the Galilean satellites and Titan may be 
accounted for by a model that delivers material to the circumplanetary disk in the 
same timescale as the feeding zone clearing timescale (< 10 5 years, see, e.g., Charnoz 
and Morbidelli 2003). In particular, the angular momentum calculation is new to 
this work 21 , as is the discussion of specific model parameters that may apply during 
satellite formation late (10 5 — 10 6 yrs) after giant planet formation. It is possible that 
inelastic collisions between planetesimals within the Hill radius of the giant planet 

form an extended (~ Rh/2) circumplanetary satellitesimal swarm, which then sus- 
21 In the present planetesimal capture framework. MEa provides a related argument in the case 

that the angular momentum budget was dominated by the gas component. Other satellite formation 

models have yet to provide an account of the angular momentum budget of the satellites of Jupiter 

and Saturn. 
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tains itself by collisionally capturing sufficiently small (< 1 km) interlopers produced 
by the fragmentation of planetesimals following giant planet formation. This process 
might result in a comparable surface density of solids in the subnebula as in the neb- 
ula (<t 2 ~ 1 g cm -2 ), and enough mass to form the Galilean satellites over several 
collisional cycles (N c ~ 10). 

This may work only if a significant fraction of the mass in the satellites originates 
from the outer regions of the planet's feeding zone, which means that the model may 
have difficulties delivering enough material to simultaneously account for the total 
masses of the regular satellites of Jupiter and Saturn. Yet, because collisions from 
Roche-lobe interlopers can pump-up the angular momentum of the circumplanetary 
disk whether or not the interaction results in significant mass capture, we show that 
a planetesimal collisional capture model may deliver sufficient angular momentum to 
account for the angular momentum budget of the regular satellites provided enough 
mass is also delivered (which may require that a significant fraction of the planetesimal 
mass be in small < 1 km fragments with low inclination orbits, as well as a solids 
enhanced solar nebula). Still, the feeding zone clearing timescale, which in this model 
controls the satellite formation timescale, may be too short to accomodate the result 
that Callisto is partially differentiated (Anderson et at. 1998; 2001). This is because 
only if the satellite grows sufficiently slowly (r acc ~ 10 5 — 10 6 years) can it radiate 
away its accretion energy in time to keep it from melting and differentiating fully 
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(Kuramoto and Matsui 1994). We have suggested that it may be possible to replenish 
the feeding zone of the giant planet with sufficiently small (1 — 100 m) planetesimal 
fragments brought in by gas drag, but we have not tested whether this is really 
possible. By this time the giant planet would presumably have opened a more or 
less clean gap in the gas disk, depending on the unknown turbulence parameter a, 
which makes this a difficult problem to tackle. Also, resonances may capture inwardly 
drifting planetesimals, which might prevent them from replenishing the feeding zone. 
However, collisions or stirring by larger particles may free them and allow them to 
continue their inward migration. It is important to note that the model we advance 
here relies on fragmentation to decrease the size of planetesimals and make them easier 
to capture, but we do not employ collisional erosion directly as a capture mechanism. 

Even if enough mass and angular momentum can be delivered, a planetesimal 
model still faces the further obstacle of how to constrain the model parameters. In 
particular, we have postponed modeling the planetesimal and satellitesimal mass 
spectrum, and have instead treated the size cutoffs for the planetesimals and satel- 
litesimals (ri and r 2 ) as adjustable parameters, meant only to indicate characteristic 
sizes containing a significant fraction of the solids disk mass. In the case of no gap 
in the planetesimal distribution, we found solutions such that r\ > r 2 ~ 1 km as 
well as T\ >> r 2 ~ 1 m. However, it is premature to conclude that this is a model 
requirement. 
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Although we have found regions of parameter space that might plausibly lead to 
a self-consistent model for regular satellite formation, until reliable ways are found 
to constrain the relevant parameters, it remains uncertain how meaningful this is. In 
particular, it is difficult to know what to choose for the size of the circumplanetary 
disk particles r 2 , where most of the satellitesimal mass resides given our assumed value 
for q 2 = 11/6. Also, our choices both for qi and q 2 are meant to be in line with earlier 
choices in Safronov et al. 1986 and other works. However, it is quite likely that further 
research will show that other choices are more appropriate in the present context. In 
particular, there may be more mass in small particles at the late times of interest than 
implied by our chosen values of q. It is unlikely, however, that this will invalidate the 
general outline of the present model given that small particles are easier to capture. 
One might argue that irregular satellites may provide an indication, but it is presently 
unknown whether the irregular satellites were collisionally captured, or bear any 
connection to the process that led to the formation of the regular satellites. On the 
other hand, we do not believe that observations of irregular satellite properties can be 
used to rule out the possibility either. Even if a connection can be made, subsequent 
collisional and dynamical evolution of satellitesimals would make it very difficult to 
extract useful constraints from the observed population of irregular satellites. 

In theory, it might be possible to model the coupled problem of the fragmenta- 
tion and growth of circumplanetary satellitesimals and heliocentric planetesimals, but 
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at present we are far from realizing this goal. This raises the issue of what physi- 
cal process sets the size for r 2 . If this size is determined by a competition between 
collisional disruption by heliocentric interlopers and Safronov-type growth in the cir- 
cumplanetary disk, then a difficulty arises for the model we advance here (at least 
for r\ > r 2 ). Namely, our system constraints lead us to a model in which collisions 
between satellitesimals and similarly-sized planetesimals are more frequent than col- 
lisions between similarly-sized satellitesimals. This is because we seek a mass balance 
(i.e., satellitesimal collisions remove a similar amount of mass as gets delivered by 
captured planetesimals), yet hypervelocity collisions with similarly-sized interlopers 
do not always result in capture (i.e., p(R) << 1). As a result, satellitesimals may not 
have sufficient time to grow to the typical sizes r 2 ~ 1 km we use to satisfy our system 
constraints (which require long satellitesimal collisional times ~ 10 5 years, such that 
N c = T acc /r co u ~ 10). However, it may be possible to circumvent this problem if we 
decrease the satellitesimal Safronov growth time by accreting "cold" circumplanetary 
debris, and increase p(R) by capturing sufficiently small heliocentric planetesimals in 
low inclination orbits and large Safronov parameter d\ » 1. It is also possible that 
the satellitesimal cutoff size is instead set by the sizes of the largest captured plan- 
etesimal fragments, or by a gravitational instability (which would set an upper-size 
cut-off of r 2 ~ 10 — 100 m, since the escape speed from particles of that size would 
be sufficient for stability), though we have not yet investigated the implications of 
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these possibilities. It seems likely that a low velocity dispersion population of small 
planetesimal fragments, as might be expected at the time of satellite formation, might 
be needed in any case. Perhaps the most promising region of parameter space is one 
in which both the satellitesimal and planetesimal populations of interest can be best 
described as collisional rubble, so that r 2 < 100 m (as hinted at by some of our solu- 
tions), but we leave investigation of that regime for future work. It should be noted 
that such a regime would likely also involve different q values than were employed in 
this paper. 

Much more work would need to be done to constrain this problem to the point 
where it can be stated whether a planetesimal collisional capture satellite forma- 
tion model is viable. The arguments presented here constitute only the outline of a 
promising model. For instance, it is clearly an oversimplification to assume as we have 
done here that collisions between retrograde and prograde satellitesimals control the 
evolution of the swarm in a timescale r coH . Also, the capture probability p(R) does 
not explicitly treat the physics of collisional fragmentation and may be unrealistic 
(we are currently considering other possible choices based on pi-scaling theory for 
the fractional mass captured following a disruptive collision). Furthermore, we have 
avoided the issue of the contribution to the angular momentum from regions of the 
disk (between ~ 0.5 — 0.7Rh] e.g., Nesvorny et al. 2003) where only retrograde orbits 
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are stable 22 and prograde objects are rapidly lost. At this time, it is most useful to 

discuss how Cassini observations can affect this model's prospects compared to that 

of MEa,b (see also Appendix C). The main observational objection that might be 

levelled against this model is that it fails to account for the composition of Saturn's 

icy satellites. In particular, the medium-sized Saturnian satellites have densities that 

would seem to preclude a solar composition ice/rock ratio (this is only marginally 

true for Titan itself), which one might not expect from a model that derives the bulk 

of the solids that make-up the regular satellites directly from heliocentric orbit. This 

issue has been noted before in a number of publications (e.g., Johnson et al. 1987; 

Lunine and Tittemore 1993; Podolak et al. 1993) largely on the basis of a comparison 

between the composition of regular satellites and objects thought to have formed in 

the outer solar system such as Pluto-Charon (see, e.g., McKinnon et al. 1997) as 

well as Triton (McKinnon and Mueller 1989). The recent Cassini flyby of Phoebe has 

yielded a density of 1.6 g cm' 3 for this retrograde, irregular satellite (Porco et al. 

2005), which indicates a rock to ice ratio of at least 50% (Johnson and Lunine 2005), 

and reinforces the demarcation between objects that formed in the outer solar nebula 

and the giant planet subnebula. 

22 One would need to argue that either the capture probability decreases outside ~ Rh/2 or that 

most of this retrograde material is not delivered to the inner disk where regular satellites form but 

rather escapes outward. 
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Mosqueira and Estrada (2005) sketch a possible resolution of this issue, involving 

a hypervelocity collision between Titan and a Triton-sized differentiated interloper, 

and the accretion of secondary icy moons out of a volatile-rich disk formed in the 

aftermath of the collision. In this scenario Iapetus may have been scattered by Titan 

to its present position, and Hyperion may represent a remnant from the collision 

which became captured into resonance and was thus prevented from crossing Titan's 

orbit and getting either scattered or accreted by it (which was presumably the fate of 

other objects forming in the region). If so, one would expect an icy composition for 

both Iapetus and Hyperion. 23 Regarding Titan itself, such an energetic impact would 

very likely have led to a fully differentiated state (Tonks and Melosh 1992), but there 

are a number of reasons to expect that Titan will turn out to be fully differentiated. 

Conversely, in the unlikely event that Cassini data indicated that Titan like Callisto 

is only partially differentiated that would argue against such an event, but in favor 

of a planetesimal collisional capture model, since it would provide evidence that the 

bulk of the material forming Titan was delivered over a long timescale as we assume 
23 In the model of MEa,b one also expects an icy composition for Iapetus and Hyperion. In this 

model satellitesimals forming in the outer disk are icy (Mosqueira and Estrada 2005) and drift 

in a long timescale due to gas drag. As they approach a satellite they might get captured into 

resonances, and collisions among them might grind them down and knock them out of resonance. 

Thus, in this model too Hyperion might represent a collisional remnant. Also, one might expect 

that Hyperion-like objects might be typical impactors in the formation of Titan and Callisto. 
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here. 

It is also important to consider the consequences of this model to the delivery 
of volatiles trapped in planetesimals to the atmosphere of Titan. This delivery of 
trapped volatiles is presumably also linked to the Galileo probe observation that 
Jupiter's heavy volatile element global abundances, excluding neon and oxygen, are 
about 3 times solar (Atreya et al. 1999; Mahaffy et al. 2000). Given that to form 
the clathrate of Ar, the model of Gautier et al. (2001a,b) requires a very low solar 
nebula temperature T ~ 36 K, which may be too cold for the location of Jupiter 
even for passive disk models, and the temperature constraint to trap it in amorphous 
ice is even lower T ~ 25 A" (Notesco and Bar- Nun 2005), the observed enhacement 
appears to argue that planetesimals trapped argon either as clathrate or in amor- 
phous ice in cold regions of the disk before migrating to warmer regions. If so, our 
planetesimal collisional capture formation model would seem to imply Ar/N 2 ~ 0.06 
Owen (2000) and that the molecular nitrogen in Titan's atmosphere was delivered as 
trapped N 2 . However, before such a conclusion can be reached one must first gain 
an understanding of the consequences of planetesimal collisional grinding on trapped 
volatiles. Thus, a planetesimal fragment capture model may be consistent with recent 
Cassini observations of Titan (Flasar et al. 2005b; Niemann et al 2005). 
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Appendix A: List of Symbols 



a Distance to the Sun. 

b Ratio of disk location R to Hill radius. 

c Gas sound speed. 

Ci Proportionality constant for rij. 

Cj Jacobi constant. 

Cd Gas drag coefficient. 

J- cs Mass delivery rate per unit volume for collisions 

between comparable-sized bodies. 
Ti Mass delivery rate per unit volume due to population i. 

G Gravitational constant. 

Hi Scale-height of population i. 

I Total collisional mass inflow rate per unit area. 

I sp Contribution to / due to satellitesimals smaller than incoming 

planetesimals. 

lib Contribution to / due to satellitesimals larger than incoming 

planetesimals. 
X Total mass collisional rate per unit area 

X L Total angular momentum collisional rate per unit area 

Jj Rate of objects in population i undergoing collisions per unit 

volume. 

L z z-component of the angular momentum. 

l z z-component of the specific angular momentum. 

£ Specific angular momentum boundary value. 

(£ z ) Mean specific angular momentum. 
m i,max Maximum cutoff in the mass distribution. 
m i,min Minimum cutoff in the mass distribution. 
M P Mass of the giant planet. 
M Mass of the Sun. 

Mi Mass delivery rate to the circumplanetary disk due to free-free 

collisions. 
M Mass accretion rate. 

M sw Steady-state mass of the circumplanetary disk. 
rii Mass distribution function. 

N c Collisional cycles. 

n p Number of times a planetesimal crosses the circumplanetary disk 

before being scattered. 
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Pi 


Orbital period. 


Pi 


Probability of capture due to collisions with population i. 


Qi 


Power-law index for n*. 


Ti,max 


Maximum cutoff in the size distribution. 


Ti,min 


Minimum cutoff in the size distribution. 


n 


Maximum cutoff in the size distribution. 


r c 


Characteristic size for which all particles could be captured into the 




circumplanetary disk. 


Rh 


Hill radius of the giant planet. 


Rp 


Radius of the giant planet. 


R 


Inner boundary of the circumplanetary disk. 


Rd 


Outer boundary of the circumplanetary disk. 


R 


Distance to the giant planet. 


Rgap 


Gap size in giant planet feeding zone. 


S 


Ratio of differential rotation in the circumplanetary disk 




to the random velocities of incoming planetesimals. 


T K 


Orbital period of the giant planet. 


T 


Characteristic drag time. 


V\ 


Mean velocity of incoming planetesimals. 


V 


Random velocity. 


V e 


Escape velocity from the giant planet. 


V K 


Kepler velocity of the giant planet. 


e 


Angular momentum deposition efficiency 


9i 


Safronov parameter. 


A G 


Mean free path of a gas molecule. 


P 


Ratio of secondary to primary. 


Pi 


Volume density of solids. 


Ps 


Planetesimal or satellitesimal density. 


Pg 


Nebular gas density. 


°i 


Solids surface density. 


T~zi 


Vertical optical depth when all particles are of size ri jinax . 


Tcoll 


Collision time in the circumplanetary disk. 


T~acc 


Planetesimal delivery or accretion timescale. 


T frag 


Fragmentation timescale for r 2;max . 




Planetesimal ejection timescale. 




Angular momentum distribution function. 




Orbital frequency. 
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Appendix B: 4 th Order Expansion of $(R,£ Z ) 

In this appendix, we derive the 4th order expansions for the mass and angu- 
lar momentum collision rates per unit area which are described by the distribution 
$(R,£ Z ), where &(R,£ Z ) d£ z is the mass collision rate per unit area of planetesimals 
with ^-component of specific angular momentum £ z in the range l z and £ z + d£ z . 
Our treatment follows that of DT93 for the case of weak gravity and high velocity 
dispersion. The mass collision rate per unit area 1{R) and the angular momentum 
collision rate per unit area Tl{R) correspond to the zeroth and first moments 

1(R) = 7^ J ®(R,t z )d£ z , 1 L (R) = ^J £MR,£ z )d£ z , (Bl) 

where integration is over all angular momenta such that \£ z \ > £q. 

In polar coordinates, the distribution funcion £ z ) is given by Eq. (43) where 
the quantity \l/ is 



1 _ v 2 cos 2 j) ^sin^-m^cos {x + tp)f 1 n p p . , 

2* = ^T + 2vKl ~ W) = ^o~B 1 + B 2 , (B2) 



where = (v /vi) 2 g(ip) and g{jp) is defined in Eq. (46), and B x and B 2 are given by 



E = V V ( h \ sin ^ cos (x + ^) B = !p 2 f b Y cos 2 (x + -0) (B3) 

1 vAs) i-ir 2 2 U i-jr ' 1 ' 
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We derive an expression for the distribution function <$>(R,£ Z ) oc = e~^ q>0 e Bl ~ B2 
by expanding the exponential in terms of B 1 and B 2 to fourth order in the quan- 
tity b/s. Using the trigonometric identity cos (x + = cos x cos — sin x sin and 
surpressing terms that are odd functions of sin and cos which will later vanish upon 
integration over ip, the expansion is 



e a-a s3 l_2/ 1 i(S)sin X + (i) 



+(*)' 



2r2(v_ 
?,J1 \ j;i 



2 A Gr) 2 - 
2^(^) 3 -3r/!(^ 

2r/i(^) 2 +r 



-2(!) 3 



(/21 - 2/22 sin \) 



(/21 - / 3 2sin 2 x)sinx , ( B4 ) 



(/41 - 2/42 sin 2 x + /43 sin \) 



where the functions / are given by 



= ^MHiHMEk , M^/,,(,) + t±ilWHM i . ( B5) 

Carrying out the integration over x where sinx = — VLiR 2 )/vR leads to 



2(JiT ; 



z2 



i?[27r (1 - |r)]V2 Pl ^ 



/»27T /» 

/ # / 

Jo J\t 



ve 



\e z -n 1 R 2 \/R 



A + A 2 [ — 

Vl 



A A I - 



Vl 



dv 



(B6) 
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where the Ai are functions of the functions f,T, and a = (£ z — Vt\R 2 ) / \f2v\R arrived 
at by collecting terms with common velocity powers in Eq. (B4). This expression 
can then be integrated over the velocity to yield Eq. (45), where the coefficients Ki 
are given by 



Ka = 1 



r/21 - J/1/21 (!) + |r 2 /4i - y/i/ 2 2 i + Okfifai (!) > 



K x = -2V2A (!) + 6v^r/!/ 21 - ^A 2 / 21 (!) 3 , X 2 = -4/i(2/ 22 - / 21 ) (!) 



f>\3 



+4 



r/i(2/ 42 - / 4 i) + f r 2 / 4 3 - ^/! 2 (/ 42 - Ai) (!) 4 , 



K, = 8V2f?(f32 - /2l) (!) 3 , ^4 = |/! 2 (/43 - 2/ 42 + /4l) (!) 4 , 

K 5 = 2gTf 22 (f) 2 - gT 2 f 42 Q) 4 , K 6 = -GgyfiFftfa (!) 3 , 



^7 = 



4^/43 + ^^/43 (!) 



"f>\4 



(B7) 



The mass collision rate per unit area is then given by integration of £ z ) over 
all \£ z \ > £q and yields 



ARM= , ~ r fTTwE^W).^^), (B8) 

7r (4 — 4fi P-, Jo [gWr /2i H 



with the .Mi given by 
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Mo = K h (ij), Mi 




M 2 = K 2 {^oW + (a e^ a o + a' Q e^<) } = K 2 {^oW + h 2 (^)] , 

M 3 = 5^ { (ag + J) e"^ - (a? + ±) e"^ 2 } = K 3 h 3 ^), 



M, = K 4 {^hoM + [ao (ag + j) e"** + aj, (a? + §) e"^] } 

M 6 = f {2hM +M)}, M 7 = f {^hoty) + 2/i 4 ty) 

(B9) 

with = (5 , / 7r ) 1 ^ 2 [«oEi(— ^ag) + (— l) n+1 a' n Ei(— pa' 2 )], where the exponential 

integral Ei(— go?) = —E^go?), and 



Equation (B8) can now be integrated numerically over the angle ip. Likewise, the 
z-component angular momentum collision rate per unit area is 



hoW = 1 



cr 



'f(^ 1 / 2 ao)+erf(( 7 1 /V ) 



T L {R,l Q ) 




x 



{IT ELo Mi{g(ip), R, £o) + ft / t 



EL)AW),i2,€ )} 



(BIO) 



with the Ci given by 
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£ 7 = f {4/i 5 (^+j 6 (V)} 

2 



(Bll) 



Note that in the limit of b/s « 1, and a = «o = 0, the integrals over ip in Eq. (B8) 
and (BIO) are identical to those in Eq. (44) and the expression for A = (£ z ) /QRp in 
Eq. (45) of DT93 (omitting the multiple-hits correction) with R = R P for the case 
of high dispersion and weak gravity. 
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Appendix C: Comparison Between the SEMM and GPPC Models 





SEMM Model 


GPPC Model 


Callisto 


Consistent with partial 


Consistent with partial 




differentiation 


differentiation 


Ganymede 


Consistent with full 


Consistent with full 




differentiation 


differentiation 


Titan 


Fully differentiated 


Partially differentiated (unless 






a large collision took place) 


Iapetus 


Formed in-situ, consistent 


Inconsistent with low density 




with low density 


Formed closer in? 


Hyperion 


Icy satellitesimal 


Icy remnant of 




fragment 


a giant collision 




resonance capture 




Low density of 


Preferential loss of 


Possible Titan collision 


inner Saturnian 


silicates 


with differentiated 


satellites 




interloper 


Compositional 


Due to temperature 


Laplace resonance 


gradient of Galilean 


gradient of subnebula 


and stochastic 


satellites 




collisional component 


Titan's eccentricity 


Possibly satellite disk 


Possible collision 




interactions 


with Triton-sized 






interloper 


Titan's isolation 


Survival by 


Break-up deep in 




gap-opening 


planetary potential 






well 


Titan's hydrocarbon 


Volatile rich 


Volatiles present 


rich atmosphere 


satellitesimals 


in planetesimals 


Size of disk from 


Specific angular 


Specific angular 


which satellites 


momentum of gas 


momentum of planetesimals 



formed 
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Satellite's 
size-distance 
correlation 



Gas drag sorting 



Radial dependence 
of isolation mass 



Separation of 
regulars and 
irregulars 



Region is cleared 
by gas drag 



An abundance of ~ 1 km 
objects yet to be 
found 
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Figure 1: Plot of the satellitesimal cutoff size r 2 (solids line) and solids surface 
density a 2 (dotted line) as a function of r\ for a case in which the solar nebula is 
ten times minimum mass of solids, a planetesimal feeding time r acc = 10 6 years, and 
0\ = 4. The extremes of N c are indicated. 

Figure 2: Plot of the optical depth t z2 (solid line) and collisional time t coU (dotted 
line) as a function of r\ for a case in which the solar nebula is assumed to have ten 
times the minimum mass of solids, a planetesimal feeding time of r acc = 10 6 years, 
and 9i = 4. The extremes of N c are indicated. 

Figure 3: Same as fig. 1, except that the solids surface density in the solar nebula 
is two times minimum mass and Q\ = 75 (characteristic of a colder planetesimal 
population). In this case, the range of values for r\ are somewhat broader, but 
significantly narrower for r 2 . 

Figure 4: Same as fig. 2, except that the solids surface density in the solar nebula 
is two times minimum mass and 6\ = 75. About N c ~ 5 collisional cycles are needed 
to complete the accretion of the satellites in the case of r\ = 10 km and r 2 = 1 km. 

Figure 5: Plot of the mass accreted over 10 6 years as a function of the planetesimal 
Safronov parameter 0\. For these curves r z2 ~ 5 x 1CT 6 , r\ = 10 km, and r 2 = 1 km. 
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Figure 6: The steady-state disk mass M sw as a function of collisional cycles N c for 
r acc = 10 6 years. The mass M sw is removed to the inner disk each r coU so this mass 
must be replenished each cycle. 

Figure 7: Mass (dotted lines) and angular momentum (solid lines) in units of the 
mass and angular momentum of the Galilean satellites (M sats ~ 4 x 10 26 g, L Z:Sats ~ 
4 x 10 43 g cm 2 s -1 ) accreted by the circumplanetary disk in r acc = 10 6 years as 
a function of the gap size R gap = y/£o/Qi. For this case N c = 10, implying a 
steady-state disk mass is M sw = 0.1M sats . The proper budget of mass and angular 
momentum to account for the Galilean satellites is achieved where the solid and the 
dotted lines first intersect at L z /L Z:Sats = M/M sats = 1. For example, for the case of 

(Ji = 4.5<7 M M, Rgap ~ 

Figure 8: Specific angular momentum delivered in units of the specific angular 
momentum of the Galilean satellites {£ z ) sats ~ 1 x 10 17 cm 2 s _1 as a function of gap 
size Rgap = v^o/^i- These curves correspond to cases such that (£ z ) / (£ z ) sats — 1 
matches the mass and angular momentum budget of the Galilean satellite system for 
R gap ~ 0.5 - 1.5R H . 

Figure 9: Specific angular momentum delivered in units of the specific angular 
momentum of the Galilean satellites (£ z ) sats ~ 1 x 10 17 cm 2 s" 1 as a function of 
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circumplanetary disk size Rd for three different values of the gap size R gap = y/io/tli, 
Rga P = Rh (solid line), R gap = 0.5R H (dotted line), and R gap = (dashed line). 

Figure 10: Mass (dotted lines) and angular momentum (solid lines) in units of 
the mass and angular momentum of the Galilean satellites (M sats ~ 4 x 10 26 g, 
L z ,sats ~ 4x 10 43 g cm 2 s -1 ) accreted by the circumplanetary disk in r acc = 10 6 years as 
a function of the gap size R gap = 5R . For this case N c = 10, implying a steady-state 
disk mass is M sw = 0.1M sats . The proper budget of mass and angular momentum 
to account for the Galilean satellites is achieved for the cases of e = 0.5 and e = 0.8 
where the solid and the dotted lines first intersect at L z /L z ^ sats = M/M sats = 1. No 
solution was found for the case e = 0.2. 

Figure 11: Specific angular momentum delivered in units of the specific angular 
momentum of the Galilean satellites (£ z ) sats ~ 1 x 10 17 cm 2 s _1 as a function of gap 
size R gap = 5i? - The curves for e = 0.5 and e = 0.8 correspond to cases such that 
(£z) I (^z) sats = 1 matches the mass and angular momentum budget of the Galilean 
satellite system for R gap ~ 0.5 — 1.5Rh- No solution was found for the case e = 0.2. 
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Figure 4: Estrada and Mosqueira 
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